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The applicability of the Fisher equation, which combines diffusion with logistic nonlinearity, 
to population dynamics of bacterial colonies is studied with the help of explicit analytic solutions 
for the spatial distribution of a stationary bacterial population under a static mask. The mask 
protects the bacteria from ultraviolet light. The solution, which is in terms of Jacobian elliptic 
functions, is used to provide a practical prescription to extract Fisher equation parameters from 
observations and to decide on the validity of the Fisher equation. 



PACS numbers: 87.17.Aa, 87.17.Ee, 87.18. Hf 



INTRODUCTION 

Evolution of bacterial colonies is a subject of obvi- 
ous medical importance and has been studied recently 
[jj], ||, [| ||, D experimentally as well as theoretically. 
Some theoretical descriptions have avoided phenomena 
such as mutation and have focused on growth, competi- 
tion for resources, and diffusion. In terms of the respec- 
tive parameters a (growth rate), b (competition param- 
eter), and D (diffusion coefficient), the basic equation 
governing the spatio-temporal dynamics of the bacterial 
population u (x, t) at a position x and time t has been 
taken to be the Fisher equation M 



du(x,t) ^d 2 u(x,t) 



dt 



dx 2 



+ au (x, t) — bu 2 (x, t) . (1) 



For simplicity, we consider throughout this paper only 
the 1-dimensional situation which is indeed appropriate 
to the observations reported. Experiments have been car- 
ried out with moving masks |3|] and observations have 
been reported about extinction transitions suggested ear- 
lier in theoretical calculations JJ] and in numerical simu- 
lations ||. Those theoretical calculations have focused on 
systems in which the growth rate a varies from location to 
location in a disordered manner, and have employed tech- 
niques based on linearization of the Fisher equation. The 
first feature has allowed the analysis to use concepts from 
Anderson localization [Q, a phenomenon well-known in 
solid state physics of quantum mechanical systems. The 
second feature has relegated the nonlinearity character 
of Fisher's equation to a secondary role. Because we sus- 
pect nonlinear features represented by —bu 2 in ([]]) to be 
of central importance to bacterial evolution, we have de- 
veloped a theoretical approach which generally retains 
the full nonlinearity of that competition term. In the 
present paper, which is the first of a series built on this 
approach of maintaining the nonlinearity in the equa- 
tion, we focus our attention on the effect of a mask on 
the spatial distribution of the stationary population of 



the bacteria. 

Consider, as in the moving mask experiments ||, an 
effectively linear petri dish in which a mask shades bac- 
teria from harmful ultraviolet light which kills them in 
regions outside the mask but allows them to grow in re- 
gions under the mask. Unlike in the moving mask exper- 
iments, however, consider that the mask does not move 
but is left stationary. Interest is in the x -dependence of 
the stationary population of the bacteria. As in previous 
considerations ||, we will assume that the growth rate 
has a positive constant value a inside the mask, and a 
negative value outside the mask. 

If we take the value of a outside the mask to be neg- 
ative infinite to reflect extremely harsh conditions (due 
to ultraviolet light) when the bacteria are not shaded 
from the light, we can take the population at the mask 
edges and outside to be identically zero. We will put 
du(x,t) /dt — in (Q ) to reflect stationarity, introduce 
a scaled position variable £ = x/\/D for simplicity, and 
begin our analysis with the ordinary differential equation 
for the stationary population u (£): 



d£ 2 



bu 2 (£) = 0. 



(2) 



Our interest is in the regions in the interior of the mask of 
width 2w, i.e., for — w < x < w, the boundary conditions 
being u \ ±w/\/D \ = 0. 

The purpose of our investigation is to give a practical 
prescription to decide on the applicability of the Fisher 
equation to bacterial evolution, and to extract the param- 
eters D, a, b from observations if the equation is found 
to be applicable. 



ELLIPTIC SOLUTIONS IN THE INTERIOR AND 
EXTRACTION OF FISHER PARAMETERS 

The solutions of @ can be written in terms of Jaco- 
bian elliptic functions as follows. It is known |8) that 



2 



the square of any of cn(£, fc), sn(£, fc) or dn(£, fc), satisfies 
an equation resembling (||). Here, we use the notation 
that k is the elliptic parameter |9) rather than the ellip- 
tic modulus which is the square of k. Thus, y — sn 2 (£, k) 
is known to satisfy 



fy 



+ 4 (1 + k 2 ) y - 6fcV 



(3) 



Comparison of (||) with (||) shows that the signs of the 
linear and quadratic coefficients are the same in the two 
equations but (||) has an extra constant term on the right 
hand side. This difference, as well as the fact that the 
bacterial system has more independent parameters than 
the single k that appears in (||) , suggests that we augment 
sn 2 (cc, k) by phase and amplitude parameters, i.e., take 
as the solution of (@) within the mask 



= asn 2 (J3£ + S,k)+J, 



(4) 



and obtain the quantities a, (3, S, 7 by differentiating (Q) 
or by other means. The suffix i represents the interior 
of the mask. Symmetry considerations, specifically the 
requirement that the maximum of Mi (£) be at £ = 0, lead 
to an evaluation of S as half the period of sn 2 . A shift 
identity allows the rewriting of (0) as 



(0 =acd 2 (/?£,*:)+ 7, 



(5) 



the cd function || being simply the ratio cn/sn. 

On differentiating (||) twice w.r.t. x, using the rela- 
tionships among the elliptic functions, and substituting 
in pj), we find: 

4/3 2 (fc 2 + l)-a + 26 7 = 0, 
Qk 2 2 ~ba = 0, 
2a^ 2 (l-fc 2 )+ 7 (a- 7 6) = 0. 

Solution of this algebraic system leads to the result that 
a and 7 are proportional to each other through a factor 
which is a function only of the elliptic parameter, 



(fc 2 + 1) + VI - k 2 + fc 4 
3k 2 



We also find explicit connections between the quantities 
a, (3 and two of the Fisher parameters of the bacterial 
system a, b, 



a 



(|f) k 2 (l-k 2 + fc*) 



-1/2 



-1/2 



(6) 



P 2 = (f)(l-fc 2 + fc4, 

This allows us to write the stationary solution as 
«i (0 = (<*/&) [ f a (k) cd 2 (Vof/3 (fc) C, fc) + / 7 (fc)] (7) 



explicitly in terms of the Fisher parameters a, 6 and three 
functions of fc alone: 

4N-1/2 



/„(*) = (3/2)A: 2 (fc' 2 + fc 4 ) 
/^(fc) = (l/2)(fc' 2 + A: 4 ) 
/ 7 (fc) = (1/2) [l - (fc 2 + 1) (fc' 2 + fc 4 ) 



-1/2' 



(8) 



Here fc' 2 = 1 - fc 2 . 

Equation (j^) provides us with the means to meet 
the primary goal of this investigation. The practical 
prescription we seek for investigating the applicability 
of the Fisher equation begins with fitting (|7|) to the 
observed stationary profile. A least-squares procedure 
yields a, 6, fc. For sensitivity purposes we use the nome 
q = cxp (-wK'/K) for fitting fuj rather than fc. The 
relation 



f a (fc) + A (fc)] 



(9) 



fc 
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(fc 



U4N1/2 



(fc' ; 



fc 



4N-1/2 



between the maximum value of the bacterial population 
u m and the extracted parameters provides a check on the 
procedure. The determination of the diffusion constant 
D follows the determination of fc. For this we can use 
the boundary condition mentioned above, that van- 
ishes at the edges of the mask: £ = ±w/yD. Equation 
(Q) leads then to an implicit expression which yields the 
diffusion constant D : 



cn' 



2 ((a/4)(l-fc 2 + fc 4 ) 1/2 w/VD,k 
((fc 2 + l) - (l-fc 2 + fc 4 ) 1/2 )(l-fc 2 ) 



(10) 



fc 2 (2-fc 2 + (l-fc 2 + fc 4 ) 1/2 ) 

Our prescription for the extraction of Fisher parame- 
ters D, a, b is, thus, complete provided we can assume 
the conditions outside of the mask to be harsh enough 
to put u at the edges to vanish. This assumption can be 
tested from the observations. The question of the very 
applicability of the Fisher equation to the bacterial sys- 
tem can be addressed by the quality of the fits of solution 
to the data. Fits of poor quality would necessitate a re- 
thinking of the quadratic nonlinearities assumed in the 
equation, indeed of the entire form of the equation. 

We illustrate our practical prescription in Fig. [l]. We 
have considered two hypothetical cases of the observed 
stationary profile of he bacterial population. One per- 
tains to a situation in which the Fisher equation is ap- 
plicable: (la); the other in which it is not: (lb). The 
'data' correspond, respectively, to stationary solutions of 
(|l|) and of the so-called Nagumo equation 111]] 



du d 2 u 



C) (Au ~ Bu 2 



(11) 
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FIG. 1: Procedure to determine applicability of Fisher equa- 
tion and/or to extract parameters from observations. Shown 
is a least-squares fitting of the analytic solution, Eq(^), of the 
Fisher equation, to numerically generated 'data' by adding 
noise to theoretical predictions in two cases. In a) the Fisher 
equation can be considered applicable while in b) it cannot. 



noise having been added in each case to simulate experi- 
ments. 

The numerically generated data are plotted with circles 
while the full line curve shows the best fit. We see that 
in Fig. la, the Fisher solution matches well the data. By 
contrast, the fitting procedure fails in lb. The intrinsic 
non-linearities in the data of lb are different from those 
characteristic of the Fisher equation (compare Eqs. (Q) 
and (pi])). Some of the data features in lb, as for exam- 
ple the change in concavity and the zero derivative at the 
borders of the mask can not be reproduced by the ana- 
lytic solution Eq (0). Thus, we have shown here how one 
would determine clearly the applicability of the Fisher 
equation to a given set of observations. 

How would one proceed if, in the light of experiment, 
the Fisher equation turns out to be inapplicable in this 
way? We suggest an additional prescription to obtain 
the form of the nonlinearity from the stationary mask 
observations. The observed stationary bacterial profile 
is Ui (x) . A numerical differentiation procedure can be 



made to produce d 2 Ui (x) jdx 1 . A plot of d 2 Ui (x) /dx 2 
versus it, (x) , the different points corresponding to dif- 
ferent values of x, would either confirm Fisher behavior 
or point to nonlinearities, such as that in the Nagumo 
equation, other than that assumed in the Fisher equa- 
tion. Fig. 2 illustrates this prescription in the context 
of the assumed observations in Fig. [l]a and Fig. [j]b. 
The 'data' were numerically differentiated in each case 
and the second spatial derivative was plotted versus u as 
shown Jl2T|. 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 



-0.2 



u"(x) 

0.0 



-0.2 



1 1 1 1 1 1 1 1 1 1 r 




0.0 



- -0.2 



,o°° ° ° 



O o o oo' 



0.0 



-0.2 



0.0 0.2 0.4 0.6 0.8 1.0 

u(x) 



FIG. 2: Procedure to extract from experiment the type of 
non-linearity in bacterial evolution. Numerically obtained 
second derivative of u is plotted against u in the two cases 
a and b of Fig. 1. 

While the quadratic nonlinearity characteristic of the 
Fisher equation is compatible with Fig 2a, the curvature 
of the data in Fig. 2b immediately points out the in- 
compatibility with the Fisher equation and suggests a 
Nagumo-like alternative. 



DEPENDENCE ON MASK SIZE 

Obviously, good experimental practice should use for 
the extraction of the Fisher parameters not a single mask 
but masks of varying sizes. It is clear that the peak value 
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of the profile, u m , will decrease as the mask size is de- 
creased (alternatively as the diffusion coefficient is in- 
creased). However, what is the precise dependence of the 
stationary profile on the size of the mask, as the size is 
varied? In answering this question, one finds that the pe- 
culiarities of the elliptic solutions produce a bifurcation 
behavior: there is a minimum mask size below which bac- 
teria cannot be supported because they diffuse into the 
harsh regions where they die. We suggest that this ef- 
fect, known in the study of phytoplankton blooms [ O j, 
be used to validate the Fisher equation in bacterial pop- 
ulation as follows. 

The dependence of the peak value of the stationary 
bacterial population on k is in (||) whereas the depen- 
dence of the mask width 2w on k is obtained by inverting 
© 



sfD 

V [((fc 2 +i)-(i-fc 2 +fc 4 ) 1/2 )(i-fc 2 )1 

\ fe 2 (2-fc 2 + (l-fc 2 +fc 4 ) 1 / 2 ) 



Cll 



1/2 



(12) 

The conjunction of @ and ( |l^ ) yield the dependence 
of the profile peak on the mask size. For a given set 
of Fisher parameters, a decrease in the mask width 2w 
from large values causes a decrease in k. This decrease is 
monotonic. The value k = is reached at a finite value 
of the width. In this limit, the elliptic function cd(/3£, k) 
becomes its trigonometric counterpart cos(/3£), and jlo| ) 
reduces to 



, 2 (|^75 



0.5. 



Thus, there is a critical size 2w r of the mask: 



2w c = 7T\ 



(13) 



(14) 



No stationary bacterial population can be supported be- 
low such a size. An excellent experimental check on the 
applicability of the Fisher equation could be the deter- 
mination of this bifurcation behavior. On the basis of 
quoted || || values D w 10~ 5 cm 2 /s, a w 10~ 4 /s, we 
obtain the critical mask size to be of the order of half a 
cm, a limit that should be observable. 

If we relax the condition that the environment outside 
the mask is harsh enough to ensure zero population of the 
bacteria, Dirichlet boundary conditions used in the pre- 
vious analysis are not appropriate. In the steady state, 
the bacterial concentration just outside the borders of 
the mask would then be different from zero as a result 
of finite diffusion. While the elliptic function solution in 
Eq. (7) (but without the Dirichlet boundary condition) 
is appropriate inside the mask, it turns out to be exceed- 
ingly difficult to find a solution outside the mask. If one 
starts out with the same (elliptic) form of the solution 
outside but with a negative but finite value of a, one gets 




2w (cm) 



FIG. 3: Reduction of the critical size of the mask as a result 
of the finiteness of a outside the mask. Shown as (a) is the 
dependence of the maximum of the profile, u m , on the width 
of the mask. For comparison we give (b), the u m dependence 
on 2w in the Dirichlet case. For case (a), the inset shows the 
actual profiles for several values of the width. 



the requirement that u (£) be negative. This is not allow- 
able since u (£) is a bacterial density which must remain 
positive. Other known solutions 



u(0 = - 



(3/2) (a/b) 
cosh 2 (^) 



(15) 



are also rejected on account of their patent negativity. It 
is possible, however, to obtain reasonable solutions jl3| if 
it is assumed that the bacterial densities outside the mask 
are so small that the quadratic term proportional to b 
may be neglected in the Fisher equation for the analysis 
in the exterior of the mask. Such an analysis leads to a 
smaller critical size relative to that in ( |l4| ) . Fig. 3 shows 
the dependence of u m on mask size for both the cases 
of (a) infinite and (b) finite (b) (negative) a outside the 
mask. The inset shows the x dependence of the solution 
for the latter case. 



REMARKS 

Our interest in the present paper being in the deter- 
mination of the applicability of the Fisher equation to 
experiments currently being conducted on the bacterial 
evolution in Petri dishes, we have displayed the explicit 
solution (7) to the Fisher equation (2) in the infinite-time 
limit when a stationary mask of given width shades the 
bacteria under it from harsh conditions outside it. Such 
stationary mask experiments we propose are easier and 
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more direct for the purposes of determination of the va- 
lidity of the Fisher equation, and for the extraction of 
the parameters of the equation. It is our suggestion that 
parameters extracted in this manner may be used subse- 
quently for the analysis of moving mask experiments |}|, 
with greater confidence in the reliability of the parameter 
values. 

We have indicated explicitly how the extraction of the 
Fisher parameters may be carried out. The numerical fit- 
ting procedure in Fig. la shows the parameters relevant 
to the hypothetical observations to be D = 10 _5 cm 2 /s 
, a = 10~ 4 /s, b = 10 -8 cm 3 /s and w = 11 cm, while 
the nome q = 0.8071 [@. The procedure does produce 
parameter values when applied to Fig. lb but the qual- 
ity of the fits is poor. Such a situation would signal 
the inapplicability of the Fisher equation. The 'data' in 
Fig. lb have been generated from the Nagumo equation 
whose intrinsic nonlinearities are incompatible with those 
of the Fisher equation as is visually clear from the best 
fits. We have shown in Fig. 2 how general manipula- 
tions of the observed data may be used to suggest the 
particular form of nonlinearity to be used in the model. 
We have also concluded that the critical size effect which 
arises directly from the solution (7) is probably within 
observable limits for bacterial evolution, the size we pre- 
dict in light of quoted parameters being of the order of 
0.5 cm. This conclusion would necessitate modification 
if the actual values of D and a are different from those 
currently believed. 

In forthcoming publications we will report our analyses 
of the spatio-temporal behavior of the bacterial popula- 
tion of relevance to time-dependent experiments. Our 
basis is the Fisher equation Jl4|] and several interesting 
alternatives such as a formalism in which diffusion is neg- 
ligible but coherent motion is present [pj| , and a formal- 
ism in which long-range competition interactions produce 
an influence function and consequently striking patterns 
||l6| in bacterial populations. 
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